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ABSTRACT 



We present three newly discovered sub-Jupiter mass planets from the SuperWASP survey: WASP-54b is a heavily bloated planet of 
mass 0.636^ jj 2 ^ Mj and radius 1.653^-™ Rj- It orbits a F9 star, evolving off the main sequence, every 3.69 days. Our MCMC fit of 
the system yields a slightly eccentric orbit (e = 0.067^^11) for WASP-54b. We investigated further the veracity of our detection of the 
eccentric orbit for WASP-54b, and we find that it could be real. However, given the brightness of WASP-54 V= 10.42 magnitudes, we 
encourage observations of a secondary eclipse to draw robust conclusions on both the orbital eccentricity and the thermal structure 
of the planet. WASP-56b and WASP-57b have masses of 0.57 L™ 34 Mj and 0.672+^ Mj , respectively; and radii of \.092 +0 ^ 
Rj for WASP-56b and 0.9 1 6^Zq oil R J for WASP-57b. They orbit main sequence stars of spectral type G6 every 4.67 and 2.84 days, 
respectively. WASP-56b and WASP-57b show no radius anomaly and a high density possibly implying a large core of heavy elements; 
possibly as high as ~50 M e in the case of WASP-57b. However, the composition of the deep interior of exoplanets remain still 
undetermined. Thus, more exoplanet discoveries such as the ones presented in this paper, are needed to understand and constrain giant 
planets' physical properties. 

Key words, planetary systems - stars: individual: (WASP-54, WASP-56, WASP-57, GSC 04980-00761) - techniques: radial velocity, 
photometry 
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1. Introduction 

To date the number of extrasolar planets for which precise mea- 
surements of masses and radii are available amounts to more 
than a hundred. Although these systems are mostly Jupiter- 
like gas giants they have revealed an extraordinary variety of 
physical and dynamical properties that have had a profound im- 
pact on our knowledge of planetary structure, formation and 
evolution and unve iled the complexity of these processes (see 
iBaraffeet al.l2010[ and references the re in). Transit surveys such 
as SuperWASP dPollacco et al.ll200q) have been extremely suc- 
cessful in providing great insight into the properties of extra - 
solar planets and their host stars (see e.g.. iBaraffe et alJ uOlO). 
Ground-based surveys excel in discovering systems with pe- 
culiar/exotic characteristics. Subtle differences in their observ- 
ing strategies can yield unexpected selection effects impacting 
the emerging distributions of planetary and stellar properties 
such as orbital periods, planetary radii and stellar metallicity 
(see e.g..lCameronll201 llfor a discussion). For example WASP- 
17b dAnderson et al.ll2010bl) is a highly inflated (# p] = 1.99/?j), 
very low density plan et in a tilted/retrograde orbit, HAT-P-32b 
dHartman et alj I20TTI) has a radius equating that of its Roche 
Lobe thus possibly losing its gaseous envelope (R p \ = 2.05^j), 
and the heavily irradiated and bloated WASP- 12b ( Hebb etaD 
2009), has a Carb on rich atmosphere (Kopparapu et al.ll2012t 
Fossati et al.ll2010l) . and is undergoing atmosphe ric evaporation 



dLlama et al.ll201 It iLecavelier Pe s Etangs 2010) losing mass to 
its host star at a rate ~10~ v M) vr' 1 dLi et al.ll2010h . On the op- 
posite side of the spectrum of planetary parameters, the highly 
dense Saturn-mass planet HD 149026b is thought to have a core 
of heavy elements with ~7 0M m, needed to expla in its small ra- 
dius (e.g JSato etal]|2005L and lCarter et alj|2009l) and the mas- 



sive WASP-18b (M pl = 10Mj, Hellie r et al 



20091) . is in an or- 
bit so close to its host star with period of ~0.94 d and ec- 
centricity e = 0.02, that it might induce significant tidal ef - 
fects probably spinning up its host star dBrown et alJ 1201 lh . 
Observations revealed that some planets are larger than ex- 
pected from standa rd coreless models (e.g.. lFortnev et al] 2007. 
Bara ffe et al.1 12008) and that the planetary radius is correlated 
with the planet equilibri um temperature and anti-correlated with 
stellar metallicity (see Guillot et al. 2006; Laughlin et al.ll201lt 
lEnoch et~al1 1201 it iFaedi et al.l 1201 lh . For these systems dif- 
ferent theoretical explanations have been proposed for exam- 
ple, tidal heat ing due to unseen comp anio ns pumping up the 
eccen tricity dBodenhei mer et afl 12001: and Bodenheimer et all 
120031) . kinetic heating due to t he breaking of atmospheric 
waves ( Guillot & Sh owman 2002), enhance d atmospheric opac- 
ity (Burro ws et alJl2007h . semi c onvection dChabrier & Ba raffe 
2007b. and finally ohmic heating dBatygin et al.ll201 ll 2010; and 
Perna et aLll2012l) . While each individual mechanism would pre- 



sumably affect all hot Jupiters to some extent, they can not ex- 
plain t he entirety of the observed radi i ( Fortne v & Nettelmannl 
l2010t lLeconteetafll2010t iPerna et alJ l2012h . More complex 
thermal evolution models are necessary to fully understand their 
cooling history. 

Recently, the Kepler satellite mission released a large num- 
ber of planet candidates (>2000) and showed that Neptune-size 
candidates and Super-Earths (> 76% of Kepler planet candi- 
dates) are common around so lar-type stars (e.g., iBorucki et al.l 
1201 U and iBatalha et"aT]|2012l) . Although these discoveries are 
fundamental for a statistically significant study of planetary pop- 
ulations and structure in the low-mass regime, the majority of 
these candidates orbit stars t hat are intrinsically f aint (V > 13.5 
for ~78% of the sample of Bor ucki et al.l 1201 Tl) compared to 



those observed from ground-based transit surveys, making exo- 
planet confirmation and characterisation extremely challenging 
if not impossible. Thus, more bright examples of transiting plan- 
ets are needed to extend the currently known parameter space in 
order to provide observation constraints to test theoretical mod- 
els of exoplanet structure, formation and evolution. Additionally, 
bright gas giant planets also allow study of their atmospheres 
via transmission and emission spectroscopy, and thus provide 
interesting candidates for future characterisation studies from 
the ground (e.g. VLT and e-ELT) and from space (e.g. PLATO, 
JWST, EChO, and FINESSE). 

Here we describe the properties of three newly discov- 
ered transiting exoplanets from the WASP survey: WASP-54b, 
WASP-56b, and WASP-57b. The paper is structured as follows: 
in § 2 we describe the observations, including the WASP dis- 
covery data and follow up photometric and spectroscopic ob- 
servations which establish the planetary nature of the transiting 
objects. In § 3 we present our results for the derived system pa- 
rameters for the three planets, as well as the individual stellar 
and planetary properties. Finally, in § 4 we discuss the implica- 
tion of these discoveries, their physical properties and how they 
add information to the currently explored mass-radius parameter 
space. 



2. Observations 

The stars 1SWASP J134149.02-000741.0 (2MASS J13414903- 
0007410) hereafter WASP-54; 1SWASP J121327.90+230320.2 
(2MASS J12132790+2303205) hereafter WASP-56; and 
1SWASP J145516.84-020327.5 (2MASS J1455 1682-0203275) 
hereafter WASP-57; have been identified in seve ral northern sky 
catalo gues which provide broad-band optical (Zacharia s et al.l 
120051) and infra-red 2MASS magnitudes dSkrutskie et alJ l2006) 
as well as proper motion information. Coordinates, broad- 
band magnitudes and proper motion of the stars are from the 
NOMAD 1.0 catalogue and are given in TableQ] 



2.1. SuperWASP observations 

The WASP North and South telescopes are located in La Palma 
(ORM - Canary Islands) and Sutherland (SAAO - South Africa), 
respectively. Each telescope consists of 8 Canon 200mm f/1.8 
focal lenses coupled to e2v 2048 x 2048 pixel CCDs, which 
yield a field of view of 7.8 x 7.8 square degrees, and a pixel 
scale of 13.7" dPollacco et al.ll2006h . 

WASP-56 (V = 11.5) is located in the northern hemi- 
sphere with Declination 6 — h23h and thus it is only observed 
by the SuperWASP-North telescope; WASP-54 and WASP-57 
(V = 10.42 and V = 13.04, respectively) are located in an 
equatorial region of sky monitored by both WASP instruments, 
however only WASP-54 has been observed simultaneously by 
both telescopes, with a significantly increased observing cover- 
age on the target. In January 2009 the SuperWASP-N telescope 
underwent a system upgrade that improved our control over the 
main sour ces of red noise, such as temperatur e-dependent focus 
changes (Barros et al. 20 THlFa"edi et alJfeOl lh . This upgrade re- 
sulted in better quality data and increased the number of planet 
detections. 

All WASP data for the three new planet-hosting stars were 
processed with the custom-built reduction pipeline described in 
Pollacco et al. (2006). The resulting light curves were analysed 
using our implementation of the Box Least-Squares fitting and 
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SysR e m de-trending alg o rithms (see ICollier Cameron etall 
2006; iKovacs et alJ 120021: iTamuz et alJ l2005h . to search for 
signatures of planetary transits. Once the candidate planets were 
flagged, a series of multi-season, multi-camera analyses were 
performed to strengthen the candidate de tection. In addition 
different de-trending algorithms (e.g., TFA. iKovacs et al] [2005) 
were used on one season and multi-season light curves to 
confirm the transit signal and the physical parameters of the 
planet candidate. These additional tests allow a more thorough 
analysis of the stellar and planetary parameters derived solely 
from the WASP data thus helping in the identification of the 
best candidates, as well as to reject possible spurious detections. 

- WASP-54 was first observed in 2008, February 19. The 
same field was observed again in 2009, 2010 and 201 1 by both 
WASP telescopes. This resulted in a total of 29938 photometric 
data points, of which 1661 are during transit. A total of 58 par- 
tial or full transits were observed with an improvement in \ 2 of 
the box-shaped model over the flat light curve of Ay 2 = -7 01, 
and signal-to-red noise value (ICollier Cameron et al.l 12006b of 
SN red = -13.02. When combined, the WASP data of WASP- 
54, showed a characteristic periodic dip with a period of P — 
3.69 days, duration T14 ~ 270 mins, and a depth ~ 11.5 mmag. 
igure Q] shows the discovery photometry of WASP-54b phase 
folded on the period above, and the binned phased light curve. 

- WASP-56 was first observed during our pilot survey in May 
2004 by SuperWASP-North. The same field was also observed in 
2006 and 2007 yielding a total of 16441 individual photometric 
observations. SuperWASP first began operating in the northern 
hemisphere in 2004, observing in white light with the spectral 
transmission defined by the optics, detectors, and atmosphere. 
During the 2004 season the phase coverage for WASP-56b was 
too sparse to yield a robust detection with only ~ 10 points 
falling during the transit phase. Later in 2006 a broad-band fil- 
ter (400 - 700 nm) was introduced and with more data available 
multi-season runs confirmed the transit detection. Over the three 
seasons a total of 14 partial or full transits were observed, yield- 
ing 300 observations in transit, with a Ax 2 = -213 improve- 
ment over the flat light curve, and SN rec i = -7.02. The combined 
WASP light curves, plotted in Figure|2l show the detected transit 
signal of period = 4.61 days, depth = ~ 13 mmag, and duration 
T\4 ~ 214 mins. 

- WASP-57b was first observed in March 2008 and sub- 
sequently in Spring 2010. A total of 30172 points were taken 
of which about 855 were during transit. About 65 full or par- 
tial transits were observed overall with a Ay 2 = —151, and 
SN re d - -6.20. Figure ^upper panel shows the combined 
WASP light curves folded on the detected orbital period of 
2.84 days. Additionally, for WASP-57b there is photometric cov- 
erage from the Qatar Exoplanet Survey (QES, Alsubai 201 lj) and 
the phase folded QES light curve is shown in Figure ^middle 
panel. In both WASP and QES light curves the transit signal was 
identified with a period ~2.84 days, duration T\n ~ 138 mins, 
and transit depth of ~ 17 mmag. 

2.2. Low S/N photometry 

Several observing facilities are available to the WASP consor- 
tium and are generally used to obtain multi-band low-resolution 
photometry to confirm the presence of the transit signal detected 
in the WASP light curves. This is particularly useful in case of 
unreliable ephemerides, and in case the transit period is such that 
follow up from a particular site is more challenging. Small-to- 
medium sized telescopes such as the remote-controlled 17 -inch 



Table 1. Photometric properties of the stars WASP-54, WASP- 
56 and WASP-57. The broad-band magnitudes and proper mo- 
tion are obtained from the NOMAD 1 .0 catalogue. 



Parameter 


WASP-54 


WASP-56 


WASP-57 


RA(J2000) 


13:41:49.02 


12:13:27.90 


14:55:16.84 


Dec(J2000) 


-00:07:41.0 


+23:03:20.2 


-02:03:27.5 


B 


10.98 ± 0.07 


12.74 + 0.28 


13.6 + 0.5 


V 


10.42 ± 0.06 


11.484 + 0.115 


13.04 ± 0.25 


R 


10.0 + 0.3 


10.7 ± 0.3 


12.7 + 0.3 


1 


9.773 ± 0.053 


11.388 + 0.087 


12.243 + 0.107 


J 


9.365 ± 0.022 


10.874 + 0.021 


11.625 + 0.024 


H 


9.135 ±0.027 


10.603 ± 0.022 


11.292 + 0.024 


K 


9.035 ± 0.023 


10.532 ± 0.019 


11.244 + 0.026 


H a (mas/yr) 


-9.8 ± 1.3 


-34.9 + 0.8 


-22.0 ± 5.4 


Us (mas/yr) 


-23.5 ± 1.2 


2.9 + 0.7 


-0.6 + 5.4 



0.04 
0.02 

0.1 I 

-0.02 
-0.04 



Orbital Phase 



0.02 I 
0.01 t- 

E « 
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-0.01 V 
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0.90 



1 .00 
Orbital Phase 



Fig. 1. Upper panel: Discovery light curve of WASP-54b phase 
folded on the ephemeris given in Table [4] Lower panel: binned 
WASP-54b light curve. Black-solid line, is t he best-fit tran- 
sit model estimated using the formalism from Mandel & A20I 
(200j). 



PIRATE telescope i n the Observatori Astronomic de Mallorca 
dHolmes et al]l201lh . together with the James Gregory 0.94 m 
telescope (JGT) at the University of St. Andrews, provide higher 
precision, higher spatial resolution photometry as compared to 
WASP, and thus have an important role as a link in the planet- 
finding chain, reducing the amount of large telescope time spent 
on false-positives. Observations of WASP-56 were obtained with 
both PIRATE and JGT, while observations of WASP-54 were 
obtained only with PIRATE. 

Multiple Markov-Chain Monte Carlo (MCMC) chains have 
been obtained for both systems to assess the significance of 
adding the PIRATE and JGT light curves to the corresponding 
dataset in determining the transit model, in particular the impact 
parameter, the transit duration, and a/R*. We conclude that for 
WASP-54 the effect is not significant, never the less, the PIRATE 
light curves were included in our final analysis presented in sec- 
tion 33.21 In the case of WASP-56 instead, because we only have 
a partial TRAPPIST light curve (see section §2.4), the full JGT 
light curve, although of lower quality, is crucial to better con- 
strain the transit ingress/egress time, impact parameter and a/R*, 
allowing us to relax the main sequence mass-radius constraint. 
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Fig. 2. Upper panel: Discovery light curve of WASP-56b phase 
folded on the ephemeris given in Table [5] Lower panel: binned 
WASP-56b light curve. Black-solid line, is t he best-fit tran- 
sit mo del estimated using the formalism from Mandel & Ago] 
(200§. 



Fig. 3. Upper panel: Discovery WASP light curve of WASP-57b 
phase folded on the ephemeris given in Table [5] Middle panel: 
QES light curve of WASP-57b. Lower panel: binned WASP light 
curve of WASP-57b. Black-solid lin e, is the best-fit transit m odel 
estimated using the formalism from lMandel & Ag ol (2002). 



2.3. Spectroscopic follow up 

WASP-54, 56 and 57 were observed during our follow up cam- 
paign in Spring 201 1 with the SOPHIE spectrograph mounted at 
the 1.93 m telescope dPerruchot etalJl2 008: Bou chv et a0 2009) 
at Observatoire de Haute-Provence (OHP), and the CORALIE 
spectrograp h mounted at the 1 .2 m Euler-Swiss tele s cope at La 
Silla Chile dBaranne et al.lll996blOueloz et al.l l2000: Pepe et al.l 
120021) . We used SOPHIE in high efficiency mode (R = 40000) 
and obtained observations with very similar signal-to-noise ratio 
(~30), in order to minimise systematic er rors (e.g., the Charg e 
Transfer Inefficiency effect of the CCD, iBouchv et all [2009). 
Wavelength calibration with a Thorium-Argon lamp was per- 
formed every ~2 hours, allowing the interpolation of the s pectra l 
drift of SOPHIE (< 3 ms" 1 per hour; see iBoisse et~al](2010h . 
Two 3" diameter optical fibers were used; the first centered on 
the target and the second on the sky to simultaneously measure 
the background to remove contamination from scattered moon- 
light. During SOPHIE observations of WASP-54, 56 and 57 the 
contribution from scattered moonlight was negligible as it was 
well shifted from the targets' radial velocities. The CORALIE 
observations of WASP-54 and WASP-57 were obtained during 
dark/grey time to minimise moonlight contamination. The data 
were processed with the SOPHIE and CORALIE standard data 
reduction pipelines, respectively. The radial velocity uncertain- 
ties were evaluated in cluding known sys tematics such as guiding 
and centering errors dBoisse et al.ll2010T) . and wavelength cali- 
bration uncertainties. All spectra were single-lined. 

For each planetary system the radial velocities were com- 
puted from a weighted cross-correlation of each spectrum 
with a numerical m ask of spectral type G2, as described in 
iBaranne et all d 19961) and lPepe et ail d2002l) . To test for possible 
stellar impostors we performed the cross-correlation with masks 



of different stellar spectral types (e.g. F0, K5 and M5). For each 
mask we obtained similar radial velocity variations, thus reject- 
ing a blended eclipsing system of stars with unequal masses as a 
possible cause of the variation. 

We present in Tables |6] [7j and [8] the spectroscopic measure- 
ments of WASP-54, 56 and 57 together with their line bisec- 
tors (V sp an)- In each Table we list the Barycentric Julian date 
(BJD), the stellar radial velocities (RVs), their uncertainties, the 
bisector span measurements, and the instrument used. In col- 
umn 6, we list the radial velocity measurements after subtract- 
ing the zero point offset to CORALIE and SOPHIE data respec- 
tively (the zero-point offsets are listed in Table |U and Table [5] 
respectively). In column 7 we also give the line bisectors af- 
ter subtracting the mean value for SOPHIE and CORALIE re- 
spectively, and finally, in column 8, the radial velocity residuals 
to the best-fit Keplerian model. The Root-Mean-Square (RMS) 
of the residuals to the best-fit Keplerian models are as follow: 
RMS = 18.9 ms- 1 for WASP-54, RMS = 19.5 ms- 1 forWASP- 
56, and RMS = 24.4ms- 1 for WASP-57. 

For all Figures presented in the paper we adopted the con- 
vention for which SOPHIE data are always represented as filled 
circles and CORALIE data are represented as open squares. 
In Figures [4] to [9] we present the RVs, V jpa „, and the residuals 
O-C diagrams for the three systems. Both CORALIE and 
SOPHIE data sets are offset with respect to the radial velocity 
zero point, ysoPHiE and ycoRALiE . respectively (see Tables|4]and 
0. We examined W span to search for asymmetries in spectral 
line profiles that could result from unresolved binarity or indeed 
stellar activity. Such effects would cause the bisector spans to 
vary in phase with radial velocity. For the three systems no 
significant correlation is observed between the radial velocity 
and the line bisector, or the bisector and the time at which 
observation were taken. This supports each signal's origin as 



4 



Faedi R: Three new exoplanets from WASP 




Fig. 4. Upper panel: Phase folded radial velocity measure- 
ments of WASP-54 obtained combining data from SOPHIE 
(filled-circles) and CORALIE (open-squares) spectrographs. 
Superimposed is the best-fit model RV curve with parameters 
from Table |4] The centre-of-mass velocity for each data set 
was subtracted from the RVs (7sophie= -3.1109 kms" 1 and 
7coralie= -3.1335kms _1 ). Lower panel: Residuals from the RV 
orbital fit plotted against time. 



Fig. 5. Upper panel: The bisector span measurements of WASP- 
54 as a function of radial velocity, values are shifted to a zero- 
mean (<V spa „>soPHiE= -29 ms" 1 , <V spm >coRALiE= 48 ms" 1 ). 
Lower panel: The bisector span measurements as a function of 
time (BJD - 245 0000.0). The bisector span shows no significant 
variation nor correlation with the RVs, suggesting that the signal 
is mainly due to Doppler shifts of the stellar lines rather than stel- 
lar profile variations due to stellar activity or a blended eclipsing 
binary. 



being planetary, rather than du e to a blended eclip sing binary 
system, or to stellar activity (see lOueloz et al. 200l]) . 

- WASP-54's follow up spectroscopy was obtained from both 
the SOPHIE and CORALIE spectrographs (see Figures [4] and 
[5]). The RMS for SOPHIE and CORALIE radial velocity resid- 
uals to the best-fit model are RMS sophie = 33.6 ms" 1 and 
RMS coralie = 8.2 m s" 1 . Typical internal errors for CORALIE 
and SOPHIE are of 10-15 ms -1 . The significantly higher RMS 
of the SOPHIE residuals is mostly due to one outlier (RV = 
134 m s" 1 ). Removing this measurement results in a RMS sophie 
= 18 ms" , which is comparable to the quoted internal error. 
This discrepant value could have resulted from observations ob- 
tained in poor weather conditions. 

- WASP-56 has radial velocity data only from SOPHIE (see 
Figures and 01. The RMS of the RV residuals to the best-fit 
model is 19.4 ms" . When removing the only discrepant RV 
value at phase 0.5 (RV = -61 ms -1 ) the overall RMS reduces 
to 12 m s" 1 , comparable to SOPHIE internal error. - Finally, for 
WASP-57 the RMS of the SOPHIE and CORALIE radial veloc- 
ity residuals to the best-fit model are RMS sophie = 22.3 ms" 1 
and /?MS coralie = 26.5 ms" 1 , respectively (see Figures [8] 
and[9]). These become 14 ms" 1 and 17.4 ms" 1 respectively for 
SOPHIE and CORALIE data sets when ignoring the two mea- 
surements with the largest errors. 



E 



1 1 1 1 1 1 1 1 1 1 1 1 1 
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T / * 
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rl 



0.0 0.2 0.4 0.6 0.8 1.0 

Orbital Phase 



Fig. 6. Upper panel: Similar to Figure |4] the phase folded ra- 
dial velocity measurements of WASP-56. The centre-of-mass 
velocity for the SOPHIE data was subtracted from the RVs 
(Tsophie= -4.6816 kms" 1 ). Lower panel: Residuals from the 
RV orbital fit plotted against time. 
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Fig. 7. Upper panel: The bisector span measurements for 
WASP-56 as a function of radial velocity, values are shifted 
to a zero-mean (<V span >soPHiE- -40 ms _1 ). Lower panel: 
The bisector span measurements as a function of time (BJD- 
245 0000.0). No correlation with radial velocity and time is ob- 
served suggesting that the Doppler signal is induced by the 
planet. 
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Fig. 9. Upper panel: Same as Figures [5] and [7] we show the bi- 
sector span measurements for WASP-57 as a function of ra- 
dial velocity, values are shifted to a zero-mean(< V j/m „>soph/e= 
-2 ms" 1 , <V ' sp an> coralie- 22 ms -1 ). Lower panel: The bisec- 
tor span measurements as a function of time (BJD - 245 0000.0). 
No correlation with radial velocity and time is observed suggest- 
ing that the Doppler signal is induced by the planet. 




Table 2. Photometry for WASP-54, WASP-56 and WASP-57 



0.2 



0.4 0.6 
Orbital Phase 



Fig. 8. Upper panel: Similar to Figures |4] and [6] for WASP- 
57. The centre-of-mass velocity for each data set was sub- 
tracted from the RVs (ysoPHiE= -23.214 kms -1 and 7coralie= 
-23.228 km s~'). Lower panel: Residuals from the RV orbital fit 
plotted against time. 



Planet 


Date 


Instrument 


Filter 


Comment 


WASP-54b 


06/04/2011 
27/02/2012 


EulerCam 
TRAPPIST 


Gunn r 
l + z 


full transit 
partial transit 


WASP-56b 


16/05/2011 
11/03/2012 


TRAPPIST 
JGT 


I + z 
R 


partial transit 
full transit 




05/05/2011 


TRAPPIST 


I + z 


partial transit 


WASP-57b 


10/06/2011 


TRAPPIST 


l + z 


full transit 




10/06/2011 


EulerCam 


Gunn r 


full transit 



2.4. follow up Multi-band Photometry 

In order to allow more accurate light curve modelling of the 
three new WASP planets and tightly constrain their param- 
eters, in-transit high-precision photometry was obtained with 
the TRAPPIST and Euler telescopes located at ESO La Silla 
Observatory in Chile. The TRAPPIST telescope and its char- 
acteris tics are described in iJehin et all d201 ll) and iGillon etafl 
d201 ll) . A detailed description of the physical charac teristics and 
instru mental details of EulerCam can be found in iLendl et all 
(I2012h . 

All photometric data presented here are available from the 
NStED database Q One full and one partial transits of WASP- 
54b have been observed by EulerCam in 2011 April 6 and 
TRAPPIST in 2012 February 26, respectively. Only a partial 
transit of WASP-56b was observed by TRAPPIST in 201 1 May 
16, and a full transit was observed by JGT in 2012 March 



http://nsted.ipac.caltech.edu 
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0.96 0.98 1.00 1.02 1.04 

Orbital phase 

Fig. 10. Euler r-band and TRAPPIST 7 + z'-band follow up 
high signal-to-noise photometry of WASP-54 during the transit 
(see Table |2]). The TRAPPIST light curve has been offset from 
zero by an arbitrary amount for clarity. The data are phase-folded 
on the ephemeris from Table|4] Superimposed (black-solid line) 
is our best-fit transit m odel estimated using the formalism from 
lMandel& Agol(l2002l) . Residuals from the fit are displayed un- 
derneath. 



11. A partial and a full transit of WASP-57b were captured 
by TRAPPIST on the nights of 2011 May 5 and June 10 re- 
spectively, while a full transit of WASP-57b was observed with 
EulerCam in 201 1 June 10. A summary of these observations is 
given in Table [2] 

We show in Figures [TUl [TTJ and fT2l the high S/N follow up 
photometry (EulerCam and TRAPPIST) for WASP-54b, WASP- 
56b and WASP-57b respectively. In each plot we show the differ- 
ential magnitude versus orbital phase, along with the residual to 
the best-fit model. The data are phase folded on the ephemerides 
derived by our analysis of each individual object (see § 13.2b . In 
Figures[T0landfT2lsome of the light curves are assigned an arbi- 
trary magnitude offset for clarity. 

2.5. TRAPPIST 7 + z'-band photometry 

TRAPPIST photometry was obtained using a readout mode of 
2x2 MHz with lxl binning, resulting in a readout time of 
6.1 s and readout noise 13.5 e~pix~', respectively. A slight defo- 
cus was applied to the telescope to optimise the observation effi- 
ciency and to minimise pixel to pixel effects. TRAPPIST uses a 
special T+z' filter that has a transmittance > 90% from 750 nm 
to beyond 1 100 nm. The positions of the stars on the chip were 
maintained to within a few pixels thanks to the 'software guid- 
ing' system that regularly derives an astrometric solution to the 
most recently acquired image a nd sends pointing corrections to 
the mount, if needed (see e.g., iGillon et aLll201 ll for more de- 
tails). A standard pre-reduction (bias, dark, flat field correction), 
was carried out and the stellar fluxes were extracted from the 




0.98 1.00 1.02 1.04 

Orbital phase 



Fig. 11. TRAPPIST 7+z'-band and JGT R-band follow up high 
signal-to-noise photometry of WASP-56 during the transit (see 
Table |2). The JGT light curve has been offset from zero by an 
arbitrary amount for clarity. The The data are phase-folded on 
the ephemeris from Table [5] Superimposed (black-solid line) is 
our best-fit tran sit model estimated using the formalism from 
iMandel & Agol ((2002). Residuals from the fit are displayed un- 
derneath. 



imag es using the IR AF/DAOPHOTEI aperture photometry soft- 
ware Stetson ( 1987). After a careful selection of reference stars 
differential photometry was then obtained. 



2.6. Euler r-band photometry 

Observations with the Euler-Swiss telescope were obtained in 
the Gunn r filter. The Euler telescope employs an absolute track- 
ing system which keeps the star on the same pixel during the 
observation, by matching the point sources in each image with 
a catalogue, and adjusting the t elescope pointing between expo- 
sures to compensate for drifts dLendl et alj feon). WASP-54b's 
observations were carried out with a 0.2 mm defocus and one- 
port readout with exposure time of 30 s. All images were cor- 
rected for bias and flat field effects and transit light curve were 
obtained by performing relative aperture photometry of the tar- 
get and optimal bright reference stars. For WASP-57b no defo- 
cus was applied, and observations were performed with four-port 
readout, and 60 s exposures. Six reference stars were used to per- 
form relative aperture photometry to obtain the final light curve. 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



7 



Faedi R: Three new exoplanets from WASP 




0.96 0.98 1.00 1.02 1.04 

Orbital phase 

Fig. 12. Euler r-band and TRAPPIST 7 + z'-band follow up 
high signal-to-noise photometry of WASP-57 during the tran- 
sit (see Table 0. The TRAPPIST light curves have been offset 
from zero by an arbitrary amount for clarity. The data are phase- 
folded on the ephemeris from Table [5] Superimposed (black- 
solid lines) a re the best-fit transit m odels estimated using the for- 
malism from Mand el & Agoll d2002f) . The residuals from each fit 
are displayed underneath the relative light curves. 



3. Results 

3.1. Stellar parameters 

For all the three systems the same stellar spectral analysis has 
been performed, co-adding individual CORALIE and SOPHIE 
spectra with a typical final S/N of ~80:1. The standard pipeline 
reduction products were used in the analy sis, and the analysi s 
was performed using the methods given in Gill on et alJ §009). 
The H a line was used to determine the effective temperature 
(T e ff). The surface gravity (logg) was determined from the Ca i 
lines at 6 122 A, 6 162 A and 6439 A along with the Na i D and 
Mg i b lines. The elemental abundances were determined from 
equivalent width measurements of several clean and unblended 
lines. A value for micro-tur bulence (ft) w as determined from 
Fe i using the method of Magain (1984). The quoted error 
estimates include that given by the uncertainties in T e g, \ogg 
and f t , as well as the scatter due to measurement and atomic data 
uncertainties. The projected stellar rotation velocity (vsin/*) 
was determined by fitting the profiles of several unblended Fe i 
lines. For each system a value for macro-turbul e nce (v mac ) was 
assumed based on the tabulation bv lBruntt et al.l d2010h . and we 
used the telluric lines around 6300A to determine the instrumen- 
tal FWHM. The values for the v mac and the instrumental FWHM 
are given in Table [3] There are no emission peaks evident in the 
Ca H+K lines in the spectra of the three planet host stars. For 
each stellar host the parameters obtained from the analysis are 
listed in Tableland discussed below: 



WASP-54: Our spectral analysis yields the following 
results: r eff = 6100 ± 100 K, logg = 4.2 ± 0.1 (cgs), and 
[Fe/H]= -0.27 ± 0.08 dex, from which we estimate a spectral 
type F9. WASP-54's ste llar mass and radi us were estimated 
using the calibration of iTorres et"al] d2010h . We find no sig- 
nificant detection of lithium in the spectrum of WASP-54, 
with an equivalent width upper limit of 0.4 mA, corresponding 
to an abundance upper limit of logA(Li) < 0.4 ± 0.08. The 
non-detection of lithium together with the low rotation rate 
implied by the vsin;'*(f ro t = 17.60 ± 4.38 d), and the lack 
of stellar activity (shown by the absence of Ca n H and K 
emission), all indicate that the star is relatively old. From the 
estimated vsin/* we derived the stellar ro tation rates, and 
we used the expected spin-down timescale dBarnesI [20071) to 
obtain a value of the stellar age through gyrochronology. We 
estimate an age of 4.4*7/!} Gyr. This value also suggest the 
system is old. Although the age is not well constrained by 
gyrochoronlogy, it is in agreement with the results obtained 
from theoretical evolutionary models discussed below, which 
imply that WASP-54 has evolved off the main sequence. 

WASP-56 and WASP-57: Both stellar hosts are of spectral 
type G6V. From our spectral analysis we obtain the following 
parameters: T eff = 5600 ± 100 K, and logg = 4.45 ± 0.1 (cgs) 
for WASP-56, T eB = 5600 ± 100 K, and logg = 4.2 + 0.1 (cgs) 
for WASP-57. As before the stellar masses and radii are esti- 
mated using the ITorres et"a"D (l2010h calibration. With a metal- 
licity of [Fe/H]= 0.12 dex WASP-56 is more metal rich than 
the sun, while our spectral synthesis results for WASP-57 show 
that it is a metal poor star ([Fe/H]= -0.25 dex). For both stars 
the quoted lithium abundances take account non-loca l thermody- 
namic equilibrium corrections (Carlsson et al.l 1994f) . The values 
for the lithium abundances if these corrections are neglected are 
as follows: logA(Li) = 1.32 and logA(Li) = 1.82 for WASP- 
56 and WASP-57, respectively. These values imply an age of 
<: 5 Gyr for the former and an age of ^ 2 Gyr for the latter 
(Sest ito & Randier] 120051) . From vsin/* we derived the stellar 
rotation period P mt = 32.58 ± 18.51 d for WASP-56, imply- 
ing a gyrochronological age (Barnes 2007) for the system of 
~ 5.5^° 6 6 Gyr. Unfortunately, the gyrochronological age can 
only provide a weak constraint on the age of WASP-56. For 
WASP-57 we obtain a rotation period of P mt = 18.20 + 6.40 d 
corresponding to an age of ~ 1 .9*M Gyr. Both the above results 
are in agreement with the stellar ages obtained from theoretical 
evolution models (see below) and suggest that WASP-56 is quite 
old, while WASP-57 is a relatively young system. 

For each system we used the stellar densities p*, measured 
directly from our Marko v-Chain Monte Carlo (MCMC) anal- 
ysis (see 33.21 and also Seager & M allen-Ornelasl 120031) . to- 
gether with the stellar temperatures and metallicity values de- 
rived from spectroscopy, in an interpolation of four different 
stellar evolutionary models. The stellar density, p*, is directly 
determined from transit light curves and as such is indepen- 
dent of the effectiv e temperature determined from the spectrum 
dHebbet all 12009). as well as of theoretical stellar models (if 
Mpi <sc M* is assumed). Four theoretical mod els were used: a) 
the Pa dova stellar models (Marig o et al.ll200 8, and Girardi et al.l 
l2010h . b) the Yonsei-Ya le (YY) models dDe marque et al.l l2004l) . 
c) the Teramo models dPietrinferni etalJ 12004 and finally d) 
the V ictoria-Regina stellar models (VRSS) dvkndenBerg et al.l 
2006). In Figures [13j [14l and [15] we plot the inverse cube root 
of the stellar density p* -1 ^ 3 = R+IM* 1 ^ (solar units) against ef- 
fective temperature, T e $, for the selected model mass tracks and 
isochrones, and for the three planet host stars respectively. For 
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Fig. 1 3. Isochrone tracks from Marig o~et alJ ([2008) and 
iGirardiet all d2010t) for WASP-54 using the metallicity [Fe/H] = 
-0.27 dex from our spectral analysis and the best-fit stellar den- 
sity 0.2 p Q . From left to right the solid lines are for isochrones of: 
1.0, 1.3, 1.6, 2.0, 2.5, 3.2, 4.0, 5.0, 6.3, 7.9, 10.0 and 12.6 Gyr. 
From left to right, dashed lines are for mass tracks of: 1.4, 1.3, 
1.2, 1.1 and 1.0 M Q . 



c\i 



c 

CO 



oo 
d 




_i i i_ 



_i i i_ 



_i i i_ 



6000 



'eff 



5500 
/K 



5000 



Fig. 14. Isochrone tracks from IDemarque et al.l d2004l) for 
WASP-56 using the metallicity [Fe/H]= 0.12 dex from our spec- 
tral analysis and the best-fit stellar density 0.88 p . From left to 
right the solid lines are for isochrones of: 1.8, 2.0, 2.5, 3.0, 4.0, 
5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0 and 14.0 Gyr. From 
left to right, dashed lines are for mass tracks of: 1.2, 1.1, 1.0, 0.9 
and 0.8 M Q . 



Table 3. Stellar parameters of WASP-54, WASP-56, and WASP- 
57 from spectroscopic analysis. 



Parameter 


WASP-54 


WASP-56 


WASP-57 


1 en 


6100 


+ 100 


5600 i 


- 100 


5600 


+ 100 


Iu fe 5 


4.2 - 


i- 0.1 


4.45 h 


- 0.1 


4.2 - 


i-O.l 


A (km s"'l 


1.4 - 


i- 0.2 


0.9 + 


0.1 


0.7 - 


i- 0.2 


v Qin / * n^m ^ 

V 3111 v I IV111 j / 


a. n - 


l()8 


1 s -+- 

1 .J 3: 


n q 


3 7- 
J. I - 


h 1 3 

c 1 .J 


L A C / AA J 


— fl 77 


3: U.Uo 




Of, 


— fl 7S 


+ o 1 n 


TNa/Hl 


— fl 30 


+ n 04 


37 + 

U.J^ 3: 


o 14. 


— fl 70 


+ 07 
:£: u.u / 


L ivl 5/ 1 1 J 


— n 7 1 


d: U.UJ 


fl 74 + 


n of, 


— fl 1 Q 


:£: u.u / 




-0.16 


±0.05 


0.31 ± 


0.07 


-0.13 


±0.08 


TCa/Hl 

Lv-ci/ j. ij 


-0.15 


±0.12 


0.09 ± 


0.12 


-0.21 


±0.11 


[Sc/H] 


-0.06 


±0.05 


0.35 ± 


0.13 


-0.08 


±0.05 


[Ti/H] 


-0.16 


±0.12 


0.18 ± 


0.06 


-0.18 


±0.07 


[Cr/H] 


-0.21 


±0.12 


0.20 ± 


0.11 






[Co/H] 






0.35 ± 


0.10 






[Ni/H] 


-0.29 


±0.08 


0.21 ± 


0.07 


-0.25 


±0.10 


logA(Li) 


< 0.4 


±0.08 


1.37 ± 


0.10 


1.87: 


t0.10 


Mass (M G ) 


1.15: 


t0.09 


1.03 ± 


0.07 


1.01 : 


t0.08 


Radius (Rq) 


1.40: 


t0.19 


0.99 ± 


0.13 


1.32: 


t0.18 


Sp. Type 


F9 


G6 


G6 


Distance (pc) 


200 


±30 


255 h 


:40 


455 


±80 



Note: Mass and radius estimate using the Torr es et aU (2010) calibra- 
tion. Spectral type estimated from using the table in lGravl «|2008). 
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Fig. 15. Isochrone tracks from IDemarque et alJ d2004l) for 
WASP-57 using the metallicity [Fe/H]= -0.25 dex from our 
spectral analysis and the best-fit stellar density 1.638 p . From 
left to right the solid lines are for isochrones of: 0.1, 0.6, 1.0, 
1.6, 2.0, 2.5, 3.0, 4.0, 5.0 and 6.0 Gyr. From left to right, dashed 
lines are for mass tracks of: 1.1, 1.0, 0.9 and 0.8 M Q . 



WASP-54 and WASP-56 the stellar properties derived from the 
four sets of stellar evolution models (Table |9]l agre e with each 
other and with those derived from the Torres et al.l d2010) cali- 
bration, within their l-cr uncertainties. For WASP-57 the best- 
fit Mi, from our MCMC analysis agrees with the values de- 
rived from theoretical stellar tracks with the exception of the 
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Teramo models. The latter give a lower stellar mass value of 
0.87 ± 0.04 M Q which is more than 1-cr away from our best- 
fit result (although within 2-cr). The stellar masses of planet 
host stars are usually derived by comparing measurable stel- 
lar properties to theoretical evolutionary models, or from em- 
pirical calibrations . Of the latter, the most widely used is the 
iTorres etaD ((2010) calibration, which is derived from eclipsing 
binary stars, and relates logg and T e s to the stellar mass and 
radius. However, while T e g can be determined with high preci- 
sion, log g is usually poorly constrained, and thus stellar masses 
derived from the spectroscopic logg can have large uncertain- 
ties and can suffer from systema tics. For example t he masses 
of 1000 single stars, derived bv IValenti et alJ ([1998) via spec- 
tral analysis, were found to be systematically 10% larger than 
those derived from theoretical isochrones. A similar discrep- 
ancy was a lso found in the anal ysis of the st ellar parameters o f 
WASP-37 dSimpson et alj|201lh. W ASP-39 dFaedi et alJl20Tlh . 
and WASP-21 dBouchv et al.ll2010l) . Additionally, different sets 
of theoretical mod els might not perfectly agree with each other 
dSouthworthll2010l) . and moreover at younger ages isochrones 
are closely packed and a small change in T e g or p* can have a 
significant effect on the derived stellar age. For each planet host 
star we show a plot with one set of stellar tracks and isochrones, 
while we give a comprehensive list of the four models' results 
in Table[9] Using the metallicity of [Fe/H]= -0.27 dex our best- 
fit ste llar properties from the Padova isochrones dMarigo et alJ 
120081 and iGirardi et alJ [2Q10h for WASP-54 yield a mass of 
I.I^q'Mq and a stellar age of 6.3^4 Gyr, in agreement with the 
gyrochronological age and a more accurate estimate. The Padova 
isochrones together with the stellar mass tracks and WASP-54 
results are shown In Figure [13] According to the stellar models, 
a late-F star with [Fe/H]= -0.27 dex, of this radius and mass 
has evolved off the zero-age main sequence and is in the shell 
hydrogen burning phase of evolution with an age of 6.3*24 Gyr. 
The best-fit stellar ages from the other sets of stellar models of 
WASP-54 also agree with our conclusion. In Figure[T3]the large 
uncertainty on the minimum stellar mass estimated from inter- 
polation of the Padova isochrones is likely due to the proximity 
to the end of the main sequence kink. The Padova evolutionary 
models were selected nevertheless, because they show clearly 
the evolved status of WASP-54. 



In Figures[14j and[T5]we show the best-fit Yonsei-Yal e stellar 
evolution models and mass tracks dDemarque et alJ2004l) for the 
planet host stars WASP-56 and WASP-57, respectively. Using 
the metallicity of [Fe/H]= 0.12 dex for WASP-56 our fit of the 
YY-isochrones gives a stellar mass of l.OljMj? M Q and a stellar 
age of 6.2^ 'j 1 Gyr. This is in agreement with the Li abundance 
measured in the spectral synthesis (see Table 01, and supports 
the conclusion that WASP-56 is indeed an old system. Using 
the metallicity of [Fe/H]= -0.25 dex derived from our spectral 
analysis of WASP-57, we interpolate the YY-models and we ob- 
tain a best-fit stellar mass of 0.89+g^ M and age of 2.6+jg Gyr. 
These results also agree with our results from spectral synthesis 
and shows that WASP-57 is a relatively young system. For each 
system the uncertainties in the derived stellar densities, temper- 
atures and metallicities were included in the error calculations 
for the stellar ages and masses, however systematic errors due to 
differences between various evolutionary models were not con- 
sidered. 



3.2. Planetary parameters 

The planetary properties were determined using a simultane- 
ous Markov-Chain Monte Carlo analysis including the WASP 
photometry, the follow up TRAPPIST and Euler photom- 
etry, together with SOPHIE and CORALIE radial velocity 
measurements (as appropriate see Table [2] and Tables [6j [7] 
and [8]l. A detailed descrip tion of the method is g iven in 
ICollier Cameron et al.1 d2007l) and lPollacco et alJ d2008l) . Our it- 
erative fitting method uses the following parameters: the epoch 
of mid transit To, the orbital period P, the fractional change of 
flux proportional to the ratio of stellar to planet surface areas 
AF = R?JR\, the transit duration 7^4, the impact parameter b, 
the radial velocity semi-amplitude K\ , the stellar effective tem- 
perature r e ff and metallicity [Fe/H], the Lagrangian elements 
Ve cos u and y[e sin co (where e is the eccentricity and co the 
longitude of periastron), and the systematic offset velocity y. 
For WASP-54 and WASP-57 we fitted the two systematic ve- 
locities ycoRALiE and ysoPHiE to allow for instrumental offsets 
between the two data sets. The sum of the x 1 for all input data 
curves with respect to the models was used as the goodness- 
of-fit statistic. For each planetary system four different sets of 
solutions were considered: with and without the main-sequence 
mass-radius constraint in the case of circular orbits and orbits 
with floating eccentricity. 

An initial MCMC solution with a linear trend in the systemic 
velocity as a free parameter, was explored for the three plane- 
tary systems, however no significant variation was foun d. For 
the treatment of the stellar limb-darkening, the models o f lClaretl 
(2000. 12001 were used in the r-band, for both WASP and Euler 
photometry, and in the z-band for TRAPPIST photometry. 

From the parameters mentioned above, we calculate the 
mass M, radius R, density p, and surface gravity log g of the star 
(which we denote with subscript *) and the planet (which we 
denote with subscript p i), as well as the equilibrium temperature 
of the planet assuming it to be a black-body (r p iA=o) and that 
energy is efficiently redistributed from the planet's day-side to 
its night-side. We also calculate the transit ingress/egress times 
T12/T34, and the orbital semi-major axis a. These calculated 
values and their 1-cr uncertainties from our MCMC analysis 
are presented in Tables g] and [5] for WASP-54, WASP-56 and 
WASP-57. The corresponding best-fitting transit light curves 
are shown in Figures [1] [2] and [3] and in Figures [TOl [TT1 and[T2l 
The best-fitting RV curves are presented in Figures |4] [6] and[8] 

- For WASP-54 the MCMC solution imposing the main 
sequence mass-radius constraint gives unrealistic values for the 
best-fit stellar temperature and metallicity, as we expected for 
an evolved star. We then relaxed the main sequence constraint 
and explored two solutions: one for a circular and one for an 
eccentric orbit. In the case of a non-circular orbit we obtain 
a best-fit value for e of Q.067*" "^. This is less th an a 3-<x 
detection, and as suggested by Lucv & Sweenevl d 19711) . Eq. 22, 
it could be spurious. From our analysis we obtain a best-fit x 2 
statistic of x 2 lIC = 24.3 for a circular orbit, and^ cc = 18.6 for 
an eccentric orbit. The circular model is parameterised by three 
parameters: K, ysoPHiE and "/coralie? while the eccentric model 
additionally constrains e cos u and e sin a>. We used the 23 RV 
measurements avail able and we performed t he Lucy & Sweeney 
F-test (Eq. 27 of iLucv & Sweenevl 119711) . to investigate the 
probability of a truly eccentric orbit for WASP-54b. We obtained 
a probability of 9% that the improvement in the fit produced by 
the best-fitting eccentri city could have ar i sen by chance if the 
orbit were real circular. Luc v & Sweenevl (j 1 97 1 ) suggest a 5% 
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probability threshold for the eccentricity to be significant. From 
our MCMC analysis we obtain a best-fit value for a> = 62*?2 de- 
gree, this differs from 90 ° or 270 ° val ues expected from a n 
eccentric fit of a truly circular orbit (see iLaughlin et all 12005). 
We decided to investigate further our chances to detect a truly 
eccentric orbit which we discuss in 33.31 Table [4] shows our 
best-fit MCMC solutions for WASP-54b for a forced circular 
orbit, and for an orbit with floating eccentricity. However, based 
on our analysis in 33.31 we adopted the eccentric solution. 

- For WASP-56b the available follow up spectroscopic 
and photometric data do not offer convincing evidence for an 
eccentric orbit. The free-eccen tricity MCMC solution y ields a 
value of e = 0.098 + 0.048. The lLucv & Sweenevidl97ll) F-test, 
indicates that there is a 42% probability that the improvement 
in the fit could have arisen by chance if the orbit were truly 
circular. With only a partial high S/N follow up light curve it 
is more difficult to precisely constrain the stellar and planetary 
parameters (e.g., the time of ingress/egress, the impact parame- 
ter b, and a/R*), however the full, although noisy, low S/N GJT 
light curve (see Figure fTTTl . allows us to better constrain the 
parameters mentioned above. Therefore we decided to relax the 
main-sequence constrain on the stellar mass and radius and we 
adopt a circular orbit. 

- For WASP-57b the follow up photometry and radial ve- 
locity data allowed us to relax the main-sequence mass-radius 
constrain and perform an MCMC analysis leaving the eccentric- 
ity as free parameter. However, our results do not show evidence 
for an eccentric orbit, and the Lucy & Sweeney test yields a 
100% probability that the orbit is circular. Moreover, we find 
that imposing the main-sequence constraint has little effect on 
the MCMC global solution. Thus, we decided to adopt no main 
sequence prior and a circular orbit. 



3.3. The Eccentricity of WASP-54b 

Here we investigate possible biases in the detection of the ec- 
centricity of WASP-54b, and we explore the possibility that the 
eccentricity arises from the radial velocity measurements alone. 
It is well known that eccentricity measurements for a planet in 
a circular o rbit can only overestimate the true zero eccentricity 
dFordll2006l) . We want to quantify whether using only the radial 
velocity measurements at hand we can find a significant differ- 
ence in the best-fit model of truly circular orbit compared to that 
of a real eccentric orbit with e - 0.067, as suggested by our 
free-floating eccentric solution. Indeed, the best-fit eccentricity 
depends on the signal-to-noise of the data, on gaps in the phase 
coverage, on the number of orbital periods cove red by the data 
set, a nd the number of observations (see e.g., Zak amska et al.l 
l201lh . 

We use the uncertainty of the CORALIE and SOPHIE radial 
velocity measurements of WASP-54, the MCMC best-fit orbital 
period, velocity semi-amplitude K, and epoch of the transit To 
as initial parameters, to compute synthetic stellar radial veloci- 
ties at each epoch of the actual WASP-54 RV data set. We gen- 
er ated synthetic radial veloci ty data using the Keplerian model 
of iMurrav & Dermottl dl999h . for the two input eccentricities, 
e — and e = 0.067. We then added Gaussian noise deviates 
to the synthetic RV at each epoch, corresponding to the original 
RV uncertainties added in quadrature with 3.5 ms~' accounting 
for stellar jitter. In this way at each observation time the simu- 
lated velocity is a random variable normally distributed around 



Table 4. System parameters of WASP-54 
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Fig. 16. Histograms of the output eccentricity distributions for 
the input e — (grey solid line), and for the input eccentricity of 
0.067 (black dashed line). 



a value v( f/J + y, with dispersion ^jo 



where y is the 

centre of mass velocity. In this manner the simulated data have 
similar properties to the real WASP-54 velocities but with the 
advantage of having a known underlying eccentricity and orbital 
properties. 
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Table 5. System parameters of WASP-56 and WASP-57 



Parameter (Unit) 



WASP-56 



Value 



WASP-57 



Value 
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o 9^:9+0.063 
J.Z.Oi_o.033 
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7"pi,A=o (K) 
3 T 14 : time between 1" and 4"' contact 



We generated 1000 synthetic data sets for each input eccen- 
tricity and found the best-fit values for e. In Figure [16] we show 
the output eccentricity distributions for the input e — (grey 
solid line), and for the input eccentricity of 0.067 (black dashed 
line). Clearly, the one-dimensional distribution of the output ec- 
centricity is highly asymmetric. Because e is always a positive 
parameter, the best-fit eccentricities are always positive values. 
We used the 1000 output best-fit values of e of the two samples 
of synthetic data sets to perform the Kolmogorov-Smirnov (KS) 
test to asses our ability to distinguish between the two underlying 
distributions. In Figure[17]we show the Cumulative Distribution 
Function (CDFs) of the 1000 mock best-fit eccentricities for the 
two cases. We show in grey the CDFs for the simulated data sets 
with underlying circular orbits, and in black the CDFs for the 
eccentric case. We calculated D, the absolute value of the max- 
imum difference between the CDFs of the two samples, and we 
used tabulated values for the KS test. We are able to reject the hy- 
pothesis that the two samples have the same underlying distribu- 
tion with a confidence of 99.999%. We then conclude that the de- 
tected eccentricity of WASP-54b could indeed be real. However, 
more radial velocity observations and photometric monitoring 
during transit and secondary eclipses are needed to better con- 
strain the orbital parameters of WASP-54b. In the following, we 
adopt the eccentric MCMC model for WASP-54b. 



4. Discussion 

We report the discovery of three new transiting extra-solar plan- 
ets from the WASP survey, WASP-54b, WASP-56b and WASP- 




0.00 0.02 0.04 0.06 0.08 



0.10 0.12 0.14 0.16 0.18 
Eccentricity 



Fig. 17. Cumulative Distribution Function for the two sets of 
best-fit eccentricities. We show in grey the CDFs for the sim- 
ulated data sets with underline circular orbits, and in black the 
CDFs for the eccentric ones. 



57b. In the following we discuss the implications of these new 
planet discoveries. 

4.1. WASP-54b 

From our best-fit eccentric model we obtain a planetary mass 
of 0.634 + °^ Mj and a radius of 1.653+°-^° Rj which yields 



_ -0.024 

a planetary density of 0.141+°^ ps . Thus, WASP-54b is 
among the least dense, most heavily bloated exoplanets and 
shares similarities with low-density pl anets such as WASP- 17 b 
(lAndersonetal.l l2010bl). WASP- 31b dAnderson et all l2010d) . 
and WASP-12b dHebbet alj|2010l) . These exoplanets have short 
orbital periods, orbit F-type host stars and theref ore are highly 
irradia ted. Using standar d corel ess models from iFortnev et all 
(120071) and IBaraffe et all (120081) . we find that WASP-54b has 
a radius more than 50% larger than the maximum planetary 
radius predicted for a slightly more massive 0.68Mj coreless 
planet, orbiting at 0.045 AU from a 5 Gyr solar-type star 
(R« P ecfed=l-105 Rj ). However, WASP-54 is an F-type star 
and therefore hotter than the Sun, implying that WASP-54b is 
more strongly irradiated. The low stellar metallicity ([Fe/H] = 
-0.27 + 0.08) of WASP-54 supports the expected low plane- 
tary core-mass thus favouring radius inflation. Different mecha- 
nisms have been proposed to explain the ob served anomalously 
large planet ary radii such as t idal heating ( Bodenheimer et al.l 
I2001L 120031) . kinetic heating dGuillot & Showman! 120021) . en- 
hanced atm ospheric opacity (|Burrows et al.l 120071) . and semi- 
convection (Chabrier & Baraffe 120071) . While each individual 
mechanism would presumably affect all hot Jupiters to some de- 
gree - for example the detected non-zero eccentricity of WASP- 
54b and the strong stellar irradiation are contributing to the 
radius inflatio n - they cannot explain the entirety of the ob - 
served radii ( IFortnev & Nettelmannl l2010l IBaraffe et al.l l2010). 
and additional mechanisms are need ed to explain t he infl ated 
radius of WASP-5 4b. More recently, iBatvgin et ail d201ll) and 
iPerna et a D (I20T0I) showed that the ohmic heating mechanism 
(dependent on the planet's magnetic field and atmospheric heavy 
element content), could provide a universal exp lanation of the 
curre ntly measured radius an omalies (see also lLaughlin et al.l 
120111) . However, according to lWu & Lithwickl d2012l) Eq.6. the 
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maximum expected radius for WASP-54b, including ohmic heat- 
ing, is 1.61Rj. This value for the radius is calculated assuming 
a system's age of 1 Gyr and that ohmic heating has acted since 
the planet's birth. Therefore, if we regard this value as an upper 
limit for the expected radius of WASP-54b at 6 Gyr, it appears 
more difficult to reconcile the observed anomalously large ra- 
dius of WASP-54b (although the value is within 1-cr) even when 
ohmi c heating is considered , similarly to the c ase of WASP- 
17b (lAnderson et al.1 l2010al). and HAT-P-32b (iHartman et all 



201 ll). as discussed by IWu & Lifhwickl d2012l) . Additionally, 
Huang & Cummingl d20 1 2h find that the efficacy of ohmic heat- 



ing is reduced at high r e jf and that it is difficult to explain the 
observed radii of many hot jupiters with ohmic heating under 
the influence of magnetic drag. The ability of ohmic heating 
in inflating planetary radii depends on how much power it can 
generate and at what depth, with deeper heating able to have 
a stronger effect on the planets e v olution C Rauscher & Menou 
120 12t iGuillot & Showman! [2002b . iHuang & Cummingl i20V% 
models predict a smaller radius for WASP-54b (see their Fig. 
12). 

However, the discrepancy between observations and the 
ohmic heating models in particu lar in the planetary low-mass 
regime (e.g. iBatvgin et alJl20lTl) . shows that more understand- 
ing of planets' internal structure, chemical composition and evo- 
lution is required to remo ve assumptions l i mitin g current the- 
oretical models. Moreover. IWu & Lifhwickl d2012l) suggest that 
ohmic heating can only suspend the cooling contraction of hot- 
Jupiters; planets that have contracted before becoming subject 
to strong irradiation, can not be re-inflated. Following this sce- 
nario, the observed planetary radii could be relics of their past 
dynamical histories. If this is true, we would expect planets mi- 
grating via planet-planet scattering and/or Kozai mechanisms 
(which can become important at l ater stages of planetary for- 
mation compare d to disc migration. lFabrvckv & Trem aine 2007, 
Nagasawa et al . 2008), to show a smaller radius anomaly and 
large misalignments. This interesting possibility can be tested 
by planets with Rossiter -McLaughlin ( RM) measurements of 
the spin-orbit alignment (Rossiter 1924; ?). We use all systems 
from the RM-encyclopedia (http://ooo.aip.de/People/rheller/) to 
estimate the degree of spin-orbit (mis)alignment. We consider 
aligned every system with \A\ < 30° (a 3-cr detection from zero 
degrees: lwinn et all 1201 Ol) . and define ?; R m = (|A| - 30°)/30° 
as the measure of the degree of (mis)alignment of each system. 
This has the advantage to show all aligned systems in the region 
-1 < rjRM ^ 0. In Figure[18]we show 77RM versus the stellar tem- 
perature T e g for planets with RM measurements. We have cal- 
culat ed the radius anomaly , % as follows (R Q b s - R e x P )/Rex P , see 
also lLaughlin et all d201 lb . The latter is expressed by a colour 
gradient. We find it difficult to identify any correlation (see also 
I Jackson et al.ll2012l and their Fig. 11), however we note that not 
many planets showing the radius anomaly have measured RM 
effects and more observations are needed before any conclusion 
can be drawn. 

We investigate the radius anomaly of WASP-54b with re- 
spect to the sample of Saturn-mass planets including the latest 
discoveries and the planets presented in this work (for an up- 
dated list see the extra-solar planet encyclopaedia^. In Figure[T9l 
we plot % against stellar metallicity {upper panel), and as a func- 
tion of stellar temperature (lower panel). WASP-54b is indicated 
with a filled fuchsia circle. WASP-54b appears to strenghten the 
correlation between planets' inflated radii and stellar tempera- 
ture, and the anti-correlation with metallicity. With an irradia- 
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Fig. 18. We plot the degree of (mis)alignment, t/rm, versus stel- 
lar effective temperature for known planets. In colour we show 
the radius anomaly % (see text for details). The black solid line 
separates aligned from mis-aligned systems following our de- 
scription. Known planetary systems with characteristics across 
the parameter space are indicated as follows: W for WASP, H 
for HatNet, C for C0R0T, K for Kepler, and others by their full 
name. Our three new discoveries are indicated by fuchsia square 
symbols. We note that no RM measurement is yet available for 
WASP-54b, WASP-56b, WASP-57b. 



tion temperature of ~ 2470 K, WASP-54 b is in the temperature 
region, identified by iPerna et alj d2012l) . with T ilT > 2000 K 
(Tjrr as defined by Hen g et al. 12012 ). in which planets are ex- 
pected to show large day-night flux contrast and possibly tem- 
perature inversion, in which case the ohmic power has its max- 
imum effect. With more gas giant planet detections we can start 
to shed some light on which mechanism might be more ef- 
ficient and in which circumstan ces. For example , in the case 
of WASP-39b dFaedi et al.ll201 ll) . and WASP- 13b dSkillen et al.1 
2009; iBarros et al.l l201 1 ). two Saturn-mass planets with similar 
density to WASP-54b, but much less irradiated (T eq = 1 1 16 K, 
and T eq = 1417 K, respectively; see Figure [T9b, ohm ic heating 
could play a less significant role (see for example IPerna et al] 
2012). However, many unknowns still remain in their model 
(e.g., internal structure, magnetic field strength, atmospheric 
composition). Thus, more gas giant planet discoveries and their 
accurate characterisation are needed to compare planetary phys- 
ical properties, in order to understand their thermal structure and 
distinguish between various theoretical models. With a mag- 
nitude of V= 10.42 WASP-54 is a bright target and thus its 
spectroscopic and photometric characterisation is readily feasi- 
ble. Given the detected non-zero eccentricity we encourage sec- 
ondary eclipse observations in the IR. These observations will 
allow precise measurement of the system's eccentricity as well 
as provide fundamental information on the thermal structure of 
the planet. However, we note that even in the case of a circular 
orbit our MCMC solution for WASP-54 yields a very similar, 
inflated, planetary radius R c ^° = 1.65 Rj, see Table[4] 

Finally, WASP-54 is an old (> 6 Gyr) F9 star which has 
evolved off the main sequence and is now in the Hydrogen shell- 
burning phase of stellar evolution (see Figure [T3l). This implies 
that recently in its life WASP-54 has increased its radius by more 
than 60%, and thus it is ascending the red giant branch (RGB). 
WASP-54b thus could be experiencing drag forces, both gravita- 
tional and tidal, which will affect its orbital radius. Two main 
factors contribute to the change of the planetary orbit: 1) the 
host-star can lose mass via stellar wind which could be accreted 
by the planet resulting in an increase of the orbital radius, and 2) 
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Fig. 19. The radius anomaly, %, for planets in the mass range 
0.2< M p \ <0.7Mj versus stellar metallicity [Fe/H] (upper panel) 
and stellar effective temperature T e ff (lower panel). Data were 
taken from the exoplanet encyclopaedia. Known planetary sys- 
tems with characteristics across the parameter space are indi- 
cated as in FigureQ~8] Our three new discoveries are indicated by 
fuchsia filled circles. 



the planet's orbital angular momentum decreases due to the tidal 
drag, leading to a decrease of the orbital radiu s. This is expressed 
by a v \ = a t ide + a mass loss (IZahn|[l977L Il989h. While the second 
term is negligible for Jupiter-like planets ([Duncan & Lissauerl 
119981) . the term due to tides can become important as it is pro- 
portional to (R+/a(r)) & , where a(r) is the decreasing orbital ra- 
dius. This could re-set the clock of WASP-54b making it appear 
younger, and maybe contributing to the planet radius inflation. 
Finally, any planet within the reach of the star's radius during 
RGB and asymptotic giant branch (AGB) phases (about 1 AU 
for a Solar-type st ar) will spiral-in and eventually merge with th e 
star or evaporate dVillaver & Lividl2007l:lLivio & Sokerill984l) . 



of 0.672 Mj and a radius of 0.916 Rj. Similarly to WASP-56b, 
we find that WASP-57b has a high density (p pl = 0.873 pj) and 
small radius. As our analysis suggests that WASP-57b may be 
relatively young (~ 2.6 Gyr), it may posses a significant core- 
mass of m ore than 50 M ffi , as derived from standard evolutionary 
models (iFortnev et al.ll2007l) . Figure [19] shows that WASP-56 
and WASP-57 systems have different physical properties, for 
example the radius of WASP-57b appears to depart from the 
observed trend with stellar metallicity, and the planet possibly 
shares more similarities with the giant planet in the HAT-P-12 
system (lHartman et al.ll2009l) . The derived radii of WASP-56b 
and WASP-57b are also consistent with more re cent planetary 
models that include ohmic heating. Models from Batvgin et al.l 
(1201 lh . lWu & Lithwickl (120121) . and lHuang & Cumming| ( l2012l) . 
all agree in that planets with lower effective temperatures have 
smaller radii, and that Jupiter-mass planets with r e ff<1400 K 
exper ience no significant radius inflation at all (see e.g. Figure 
6 by Batvgin e t al.ll20lH) . However, despite all available infor- 
mation, we are still far from knowing the composition of the 
deep interior of exoplanets. For example the very existence 
of planetary cores, their masses, as well as the amount and 
distribution of heavy elements in the plane ts' core or in their 
envelo pes, remain undetermined. Recently, Wilson & M ilitzerl 
(1201 2|) suggested that planetary cores, mostly composed of 
rock and ices, can be eroded and/or dissolve (depending on 
their mass) into the metallic H/He layers abo ve, and thus be 
redist ributed in the planetary envelope (see also lUmemoto et all 
2006). This can have significant implications for giant planets' 
thermal evolution, their radius contraction, and overall structure. 

In conclusion, it is clear that continued exoplanet discoveries 
are needed to provide stronger constraints on theoretical models 
of close in giant planets and hence their physical properties. 
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4.2. WASP-56b and WASP-57b 

Our modelling of the WASP-56 system yields a planet mass of 
°- 57 l-o 035 Mj and radius of 1.092+°;°| Rj which in turn give 
a planet density of 0.438+°°^ pj Hence WASP-56b belongs 
to the cla ss of Saturn-mass pl anets and does not show a radius 
anomaly (Lau ghlin et alJl201ll) . Figure [191 shows that the radius 
of WASP-56b is not inflated. With a metallicity of [Fe/H] = 
+0.12 dex WASP-56 is more metal rich than the Sun, and with 
an age of ~6.3 Gyr s tand ard planetary evolut ionary models from 
IFortnev et al] f2007) and lBaraffe etall (120081) . show that WASP- 
56b has a core of approximately > 10 M ffi of heavy material. 
Moreover, WASP-56b orbits a main sequence G6 star, thus it is 
subject to less stellar irradiation than WASP-54b. Our best-fit 
MCMC model for the WASP-57 system yields a planetary mass 
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Table 6. Radial velocity and line bisector span measurements of WASP-54. 
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Table 7. Radial velocity and line bisector span measurements of WASP-56. 
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Table 8. Radial velocity and line bisector span measurements of WASP-57. 
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Table 9. Theoretical evolutionary models for WASP-54, WASP-56 and WASP-57. 
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